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We study current correlations in a T-junction composed of a grounded topological superconductor 
and of two normal-metal leads which are biased at a voltage V. We show that the existence of an 
isolated Majorana zero mode in the junction dictates a universal behavior for the cross correlation 
of the currents through the two normal-metal leads of the junction. The cross correlation is negative 
and approaches zero at high bias voltages as —1/E. This behavior is robust in the presence of disor¬ 
der and multiple transverse channels, and persists at finite temperatures. In contrast, an accidental 
low-energy Andreev bound state gives rise to nonuniversal behavior of the cross correlation. We 
employ numerical transport simulations to corroborate our conclusions. 

PACS numbers: 71.10.Pm, 74.45.-l-c, 74.78.Na, 73.50.Td 


I. INTRODUCTION 


Majorana bound states (MBS) in condensed matter 
physics are zero-energy modes which are bound to the 
boundaries of an otherwise gapped topological supercon¬ 
ductor (TSC). Such an MBS is described by a self-adjoint 
operator and is protected against acquiring a finite en¬ 
ergy. These properties are responsible for much of the 
great interest in MBSs^’^. 

Several theoretical proposals have been put for¬ 
ward for realizing topological superconductivity in con¬ 
densed matter systems^^^®. Promising platforms include 
proximity-coupled semiconductor nanowires®’® and fer¬ 
romagnetic atomic chains^®^^^, where recent transport 
measurements have provided compelling evidences for 
MBS formation^®^®®. 

Much emphasis has been put on investigating the dif¬ 
ferential conductance through a normal lead coupled to 
an MBS®®^®®. At low enough temperatures the differ¬ 
ential conductance spectrum shows a peak at zero bias 
voltage which is quantized to 2e®//i. The observation of 
such conductance quantization has proved to be difficult, 
because it requires the temperature to be much lower 
than the width of the peak. 

Alternatively, one can seek for signatures of an MBS 
in current correlations. Various aspects of current noise 
in topological superconducting systems have been stud- 

ied26. 30-34 

we consider a setup composed of multi¬ 
ple leads coupled to an MBS, which we term a “Majorana 
beam splitter” (Fig. 1), and study the cross correlations 
of the currents in the leads. In a recent work®® we have 
examined the cross correlation between currents of oppo¬ 
site spin emitted from an MBS, showing that it is neg¬ 
ative in sign and approaches zero at high bias voltage. 
In the present work we show that this result holds much 
more generally: The cross correlation of any two channels 
in the beam splitter has the same universal characteris¬ 
tics, i.e., it is negative and approaches zero at voltages 
larger than the width of the Majorana resonance, inde¬ 
pendently of whether the different channels are spin re¬ 
solved or not. An immediate experimental consequence 
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FIG. 1: (a) The proposed experimental setup is a T-junction 
between a topological superconductor (TSC) and two metallic 
leads. Here the TSC is realized by a semiconductor nanowire, 
proximity coupled to a conventional s-wave superconductor 
under an applied magnetic field, (b) We model the TSC by 
a spinless p-wave superconductor. It is coupled to the leads 
through a normal-metal section N, whose length dN is taken 
to zero. Scattering at the NP interface is described by the 
reflection matrix tnp [see Eq. (3)], while scattering at the 
T-junction is described by the matrix Sj [see Eq. (4)]. 


is that this effect can be observed in a much less challeng¬ 
ing setup, which does not require spin filters to resolve 
the current into its spin components. 

The rest of this paper is organized as follows. In Sec. II 
we describe the setup under study and state our main re¬ 
sults. In Sec. Ill we employ a simple model for the Ma¬ 
jorana beam splitter, and calculate the current cross cor¬ 
relation using a scattering-matrix approach. In Sec. IV 
we corroborate our conclusions in a numerical simulation 
of a microscopic model, comprising a proximity-coupled 
semiconductor wire. In Sec. V we present a semiclassical 
picture of transport, and use it to rederive our results 
in the high-voltage limit. This is done in a way which 
relates the result of this paper to the nonlocal nature of 
MBSs. Finally, we conclude in Sec. VI. 


II. SETUP AND MAIN RESULT 


We consider a T-junction between a topological super¬ 
conductor (TSC) and two normal-metal leads as depicted 
in Fig. 1(a). We study the low-frequency cross correla- 
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tion of the currents through the two arms of the junction, 
namely 


^’rl = / dt{Sin{0)Sii,{t)), ( 1 ) 

J —oo 

where <5/,, = Iri — (/?,), and /,,=r,l are the current op¬ 
erators in the right and left arm of the junction respec¬ 
tively^®. The brackets stand for thermal quantum aver¬ 
aging. We denote the width of the resonance due to the 
MBS by r, and the excitation gap by A voltage 

V is applied between the superconductor and the leads. 
Below we show that in the regime eV < A, Prl has a 
simple, universal behavior, given by Eq. (9). In particu¬ 
lar, Prl is negative, and approaches zero when eV 3> T. 
For eV > A the behavior is nonuniversal. 

This effect survives, to a large extent, at finite tem¬ 
peratures. As long as the temperature T is smaller than 
y, Prl is only weakly temperature dependent, even if 
T > r. This is in contrast to the zero-bias peak in the 
differential conductance spectrum which is only quan¬ 
tized to 2e^//i for T ^ T. 

Unlike studies which have focused on the cross corre¬ 
lation between currents through two MBSs at the two 
ends of a here the effect is due to a single 

MBS. In Ref. [30,33,38-40] it was crucial that the MBSs 
at the two ends of the TSC were coupled^^. Here, on 
the other hand, the effect is most pronounced when the 
two MBSs are spatially separated such that only a single 
MBS is coupled to the leads. 


III. SCATTERING MATRIX APPROACH 


The proposed experimental setup is described in 
Fig. 1(a). A semiconductor nanowire is proximitized to 
a grounded s-wave superconductor. When a sufficiently 
strong magnetic field is applied, the wire enters a topo¬ 
logical phase®’®, giving rise to an MBS at each end. One 
of the wire’s ends is coupled to two metallic leads, both 
biased at a voltage V. 

To calculate the currents through the leads and their 
cross correlation we use the Landauer-Biittiker formal¬ 
ism in which transport properties are obtained from the 
scattering matrix, describing both normal and Andreev 
scattering. We are interested in bias voltages smaller 
than the gap, A®^. An electron incident from one of the 
normal leads is therefore necessarily reflected from the 
middle (superconducting) leg. It can be reflected to the 
right or the left lead, either as an electron or as a hole. 
Since there is no transmission into the superconductor, 
scattering is described solely by a reflection matrix. 

Each normal lead contains 2M transverse channels, in¬ 
cluding both spin species. The overall reflection matrix 
which we wish to obtain reads 


^ tot — 
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where each block is a 4M x 4M matrix. The matrix 
element r“^, where a,f5 = {e, h}, is the amplitude for 
a particle of type /3 coming from the channel j to be 
reflected into the channel * as a particle of type a. Here, 
i = 1,, 2M enumerates the channels in the right lead 
while i = 2M -|-1,..., AM enumerates the channels in the 
left lead. 

We model the TSC as a spinless p-wave superconductor 
which is a valid description close to the Fermi energy"^®’^®. 
It is convenient to insert a (spinless) normal-metal sec¬ 
tion between the TSC and the junction. In this way, we 
separate the scattering in the T-junction itself from the 
scattering at the normal-p-wave interface (cf. Fig. lb). 
The length of the normal-metal section is then taken 
to zero. 

Andreev reflection at the normal-p-wave superconduc¬ 
tor interface is described by 

= ( “(.) "o'’) • P) 

where a(e) = exp [—f arccos(e/A)] is the Andreev reflec¬ 
tion amplitude for jej < with e being the energy 

as measured from the Fermi level. The information about 
the topological nature of the system is encoded in tnp (e). 
The relative minus sign between the off-diagonal elements 
of rNp(£) signals that the pairing potential of the super¬ 
conductor has a p-wave symmetry. Moreover, the non¬ 
trivial topological invariant^®’^^ Q = det[rNp(0)] = — 1 
dictates the existence of an MBS at each end of the su¬ 
perconductor. 

Scattering at the T-junction (which connects the added 
normal section to the two leads) is described by 


5j = 



; ^ 


e — 



( 4 ) 


where Se describes scattering of electrons and S'* de¬ 
scribes scattering of holes. Here, r is a 4M x AM matrix 
describing the reflection of electrons coming from the left 
and right leads (each having 2M transverse channels), r' 
is a reflection amplitude for electrons coming from the 
middle leg (having a single channel), t is a 1 x 4M trans¬ 
mission matrix of electrons from the right and left leads 
into the middle leg, and t' is a 4M x I transmission 
matrix of electrons from the middle leg into the right 
and left leads. The matrix Sg is assumed to be energy- 
independent in the relevant energy range, but is other¬ 
wise a completely general unitary matrix. 

We can now concatenate S'j with tnp to obtain the 
overall reflection matrix rtot of Eq. (2). The block r®® is 
obtained by summing the contributions from all the var¬ 
ious trajectories in which an electron is reflected back as 
an electron, while the block is obtained by summing 
those trajectories in which an electron is reflected as a 
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hole. This yields 

r®® = r + t'{—a)r'*at + t'{—a)r'*ar'{—a)r'*at + ... 

_ a{£ft'r'*t 

-^~l + |r'|2a(e)2-’ 

(5a) 

^he _ _j_ t'*ar'{—a)r'*at + ... 

_ a{e)t'*t (5b) 

1 + |r'pa(e)^ ’ 


The two other blocks are given by r®^(e) = [r^®(—e)]* 
and r^^(e) = [r®®(—e)]* in compliance with particle-hole 
symmetry^®. 

Given the blocks of the reflection matrix, the sum of 
currents in the leads and their cross correlation are ob¬ 
tained by"^® 


_ nOO 
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where I = (/r) -t- (/l) is the total current in the leads, 
and with/e(£) = l-/,i(-£) = l/{l-hexp[(£-eG)/fcBr]} 
being the distribution of incoming electrons in the leads. 
Here, the index i = 1,..., 2M runs only over the channels 
of the right lead, while the index j = 2AI+1, ..., 4M runs 
only over those of the left lead. We use a convention in 
which sgn(a) = 1 for a = e and sgn(a) = — 1 for a = h. 
At zero temperature Eq. (6) reduces to®° 

Op 

1=^ d£Tr(r^"r'^"t)^ 

h Jo 

= f E rden^ie), (7) 
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where 7^“^ = 

Let us introduce the parameter D = J2i=i repre¬ 
senting total normal transmission from the two leads into 
the middle leg of the T-junction. Inserting Eq. (5) into 
Eq. (7) and using the unitarity of Se, we first obtain the 
differential conductance 

dl 2e^ 

dV h (eE)2+r2’ ^ ’ 


where T = AI?/2\/l — D. As expected dl/dV has a 
peak at E = 0 which is quantized to 2e^//i. Similarly, we 
obtain for the cross correlation^® 


2p2 pV 

Prl{V) = ^ 


(9) 


where T^ = |^iP/2Vl - D (note that T = Tr -|- 

Tr). The cross correlation Prl is negative for all V and 
approaches zero as — 1/E for eE ^ T. This result is valid 
for eE < A. It is valid even in the presence of strong 
disorder in the junction region, as we did not assume a 
particular form of Sp. Moreover, it does not depend on 
a specific realization of the TSC hosting the MBS. 

The low-voltage behavior of the result in Eq. (9) can 
be understood from simple considerations based on the 
properties of MBSs. For eE <C T and at zero temperature 
the conductance through the MBS is quantized to 2e^/h, 
resulting in an overall noiseless current®^. Upon splitting 
the current into the two parts /r and Jl, the total noise P 
is related to the cross correlation via P — Pr+Pl+‘2Prl, 
where Pr and Pr are the current noises through the right 
and left leads, respectively®®. Since P —>■ 0 at low volt¬ 
age, while Pr and Pr are non-negative by definition, one 
must have Prr < 0. More specifically, at zero voltage 
the total noise obeys®® dP/dV\v=o = 0. In addition, 
since (for zero temperature) Pr(0) = Pr(0) = 0, one has 
dPB./dV\v=o,dPh/dV\v=o > 0. It therefore follows that 
dPRL/dV\v=o < 0- The cross correlation Prr is thus 
negative at low voltage. 

The high-voltage limit of Eq. (9) can be derived in a 
semiclassical picture of transport, based on the nonlocal 
nature of MBSs. In particular, the analysis relies on the 
fact that no local probe can determine the occupation of 
the MBS. This is explained in Sec. V below. 


IV. NUMERICAL ANALYSIS 

We now turn to illustrate the results of the previous 
section using numerical simulations. We consider the sys¬ 
tem depicted in Fig. 1(a). A semiconductor nanowire of 
dimensions S> Wy » lEz is proximity coupled to a 
conventional s-wave superconductor and is placed in an 
external magnetic field. 

The Bogoliubov de-Gennes Hamiltonian describing the 
nanowire is given in Nambu representation, 4'f(x) = 
by 

U + V(x, j/)]t^ -h iXRia^dj; - a'^dy)T^ 

2m, (10) 

-^Pa" + Ai„d(a;)r", 

where m, is the effective mass of the electron, E(cc, y) 
includes both the chemical potential and a disordered 
potential, Ar is the Rashba spin-orbit coupling strength, 
B is the magnetic field directed along the wire, y-B is 
the Bohr magneton, g is the Lande g-factor, Ai„d(a;) = 
lS.o9(x — xo) is the proximity-induced pair potential, and 
(T and T are vectors of Pauli matrices in spin and particle- 
hole space, respectively. Since we take lEz to be much 
smaller than the magnetic length, we can ignore the or¬ 
bital effect of the magnetic field. 


(eE)2 + r2 ’ 
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FIG. 2: (a) Zero-frequency cross correlations Prl [defined in 
Eq. (1)] of the currents through the left and right leads as a 
function of bias voltage V at various temperatures. Prl is 
negative for all V and approaches zero at voltages which are 
larger than both the resonance width and the temperature, 
(b) Total differential conductance, dl/dV, where I = /r-|-/l. 
At zero temperature dl/dV exhibits a zero-bias conductance 
peak quantized to 2e^/h,®^. A nonzero temperature widens 
the peak and reduces its height to a nonuniversal value. The 
inset shows the zero-temperature local density of states at 
zero energy in the wire in the absence of coupling to the leads 
in arbitrary units. The section of the wire not covered by the 
superconductor is a; G [0,2;o], as depicted in Fig. 1(a). 


We approximate the continuum model of Eq. (10) by 
a tight-binding Hamiltonian 

r s,s' 

- li^tbSss' + iu{(T,s' X d) • -b H.C.]} 

d—x,y 


r-x^XQ 


( 11 ) 


where r runs over the sites of an by Ny square lattice 
with spacing oq. Here ttb = l/2weaQ, u = Ai^/2ao, 14 = 
+ + fj. is the chemical potential, and 1^^'® is a 

Gaussian-distributed disorder potential with zero average 
and correlations 14 ^‘® 1 P)“ = Uj;gi5rr'. 

We express H in a, first quantized form as a AN^Ny x 
ANxNy matrix "Htb, from which one extracts the re¬ 
tarded Green function 

G^(£) = (e--HTB + i7rWWl)”\ (12) 

and subsequently the reflection matrix®^’^^ 

rtot(e) = 1 - 2mW^G^{e)W. (13) 


Here, W is a matrix describing the coupling of the eigen- 
modes in the leads to the end of the nanowire as depicted 
in Fig. 1(a) and specified in Appendix B. The metallic 
leads are described in the wide band limit by an energy 
independent W. With the help of Eqs. (2) and (6) we 
then obtain the currents through the leads and their cross 
correlation 


In the present work we use parameters consistent with 
an InAs nanowire, namely Ego = me-\]fl2 = 75/reV, Igo = 
— 130nm, and g = 20^°. The induced pair 
potential is taken to be Ag = 150^eV. The length of 
the wire is = 2/im, with the section not covered by 
the superconductor being xq = 200nm in length, and the 
width of the wire is Wy = 130nm. 

In Fig. 2 we present the cross correlation Prl(E) and 
the differential conductance dl/dV at various temper¬ 
atures for p, = 0 and B = 520mT. For these values 
of g and B the system is in the topological phase®’®’®^. 
Prl is negative and approaches zero at high voltages, in 
agreement with the analytic expression of Eq. (9). In¬ 
terestingly, this behavior persists even at nonzero tem¬ 
peratures. The main effect of temperature is to increase 
the voltage above which Prl starts approaching zero. 
Since the gap in the system is about 100/reV, the effect 
can be seen even at the relatively high temperature of 
T = lOOmK, a temperature for which the zero-bias con¬ 
ductance peak is much lower than 2e^/h. 

Next, we study the effect of disorder on Prl. Fig¬ 
ure 3(a) presents Prl for 10 different realizations of ran¬ 
dom disorder with Udis = 75/xeV. As expected, the be¬ 
havior of Prl does not change significantly. We can com¬ 
pare this to the case of an ordinary Andreev state which 
is tuned to zero energy. The end of the wire which is not 
covered by a superconductor (a; < xg in Fig. 1(a)) hosts 
Andreev bound states which are coupled to the leads. 
For each realization of disorder, we tune the magnetic 
field to bring one of them to zero energy®®, and calculate 
Prl. In all the realizations, the resulting tuned magnetic 
field was below the critical field B^ = 260mT, i.e., the 
system is in the trivial phase. As shown in Fig. 3(b), the 
behavior of Prl is nonuniversal and varies significantly 
from one realization of disorder to another. Importantly, 
in all cases Prl is positive at large V. 

In our simulations we have chosen the length of the 
wire Lx = 2/im to be sufficiently bigger than the local¬ 
ization length of the Majorana wave function (which here 
is about ^ ^ 300nm), so that the leads are only coupled 
to a single MBS. If ^ becomes of the order of L^, say by 
increasing the magnetic field B, then the leads become 
coupled also to the MBS at the other end of the wire. 
At this point it is as if the leads are coupled to a single 
ABS. Increasing the magnetic field therefore induces a 
crossover between the MBS case and the ABS case, in 
exactly the same way which was described and analyzed 
in Ref. [35]. 

It is interesting to examine the case when more than 
a single transverse channel is occupied in the wire. For 
weak pairing®^, the system is in the topological phase 
whenever an odd number of channels is occupied. Fig¬ 
ure 4 presents Prl and dl/dV for various values of /i, 
each corresponding to a different odd number of occu¬ 
pied channels between I and 7. When more than a single 
channel is occupied we can have subgap Andreev bound 
states which coexist with the MBS. One such state can be 
seen in Fig. 4(b) as a peak at H ~ 80/xeV. It is only below 
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FIG. 3: Current cross correlation Prl vs. bias voltage V at 
— 0 and T — 0 for different realization of short-range Gans- 
sian disorder, (a) B = 52QmT > Be, the system is in the 
topological phase with a zero-energy Majorana bound state 
(MBS) at each end of the wire. The universal behavior of 
Prl(F), (being negative and approaching zero at high volt¬ 
age) is not affected by the presence of disorder, (b) For each 
realization of disorder the magnetic field is tuned to have an 
Andreev bound state (ABS) with zero energy at the end of 
the wire, while keeping the system in the topologically trivial 
phase, B = 170 —200mT < Be (see the text for more details). 
The behavior of Prl(I^) varies significantly for different re¬ 
alizations of disorder. In all cases Prl > 0 for large V in 
contrast to the topological case where it goes to zero. 



(a) VIp^V] (b) VW] 


FIG. 4: (a) Cross correlation and (b) differential conductance 
at various chemical potentials fj., corresponding to a different 
odd number of occupied transverse channels. The calculations 
are performed at T = 0, Udia = 0, and B = 520mT. The 
addition of occupied channels introduces extra subgap states 
which coexist with the Majorana bound state. These appears 
as peaks in the differential conductance spectra at finite V 
[see (b) at 14 ~ SO/reV]. Above this voltage the behavior of 
Prl is no longer universal. 


this voltage that the behavior of PRL(f^) remains qual¬ 
itatively the same as in the single channel case. In this 
respect, the existence of subgap states reduces the effec¬ 
tive energy gap below which Prl(I^) exhibits its universal 
features. Another effect of introducing higher transverse 
channels is the stronger coupling of the middle leg of the 
T-junction to the two leads®®. 


V. SEMICLASSICAL PICTURE 


The behavior of the current cross-correlation, as given 
in Eq. (9), at high voltages can be derived based on sim¬ 


ple semiclassical considerations. We reconsider the setup 
shown in Fig. 1(a), and examine the limit eV ^ F, where 
r is the width of the zero-energy resonance (which can 
be either an MBS or an ABS). 

In this limit, the transport of current from the super¬ 
conductor to the leads can be described in terms of a 
sequence of tunneling events. In each tunneling event, a 
Cooper pair in the superconductor dissociates; one elec¬ 
tron is emitted into the right or left lead, and the other is 
absorbed into the zero mode localized at the edge of the 
superconductor. In the presence of such a zero mode, 
the many-body ground state of the superconductor is 
doubly degenerate. We denote the two ground states by 
|0) and |1), corresponding to an even and odd number of 
electrons in the superconducting wire, respectively. Each 
time an electron is emitted into the leads, the supercon¬ 
ductor flips its state from |0) to |1) or vice versa. 

Let us denote by T^/h and T*l/h the probability per 
unit time to emit an electron into the right or left lead, 
respectively, given that the superconductor is in state |0). 
Similarly, E^l/Zi are the corresponding rates when the 
system is in the |1) state. 

After a time r, there are IVr and Al electrons emit¬ 
ted to the right and left leads respectively. The average 
currents in the leads are given by 


T T 

and the current cross correlation is given by 

Prl = lim - [ dh [ dt2{6lR{ti)6Ii,{t2)) 
r^oo r Jo Jo 

=-((ArAl) - (AR)(iVL)). 

r 


(14) 


(15) 


In the case of a Majorana zero mode, all the local prop¬ 
erties of the states |0) and |1) are identical. This is usu¬ 
ally stated as the fact that one cannot make a local mea¬ 
surement which would reveal in which of the two ground 
states the system is in. In particular, this implies that 
= Tp = Tr and Tj* = T^ = Tl. Let us divide the 
time r into short time intervals At ~ At is the min- 
imal time between consecutive emission events (set by 
the minimal temporal width of an electron wave packet 
whose energy spread is ^ eV). At each time step At, ei¬ 
ther an electron is emitted to the right lead, an electron is 
emitted to the left lead, or no electron is emitted at all. 
The transport process is thus described by a trinomial 
distribution. The probabilities of being emitted to the 
right and left lead are pr = TjiAt/h and Pl = T RAt/h, 
respectively, and there are overall N = rjAt time steps. 
One thus obtains®® 


{Nr) = NpR = TRT/h, 

(Nr) = NpR = Yrt/H, 

(NrNr) - (Nr) (Nr) = -NprPr =-- 
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Finally, inserting Eq. (16) into Eqs. (14) and (15) one has 


and 


= ; (/L) = ^fL, (17) 


-Prl 


e fRFL 
h V 


(18) 


Prl is negative and approaches zero as —1/E. We have 
therefore reproduced the high-voltage limit of Eq. (9). 

Unlike the case of an MBS, for an ABS the probabilities 
to emit an electron to the right or the left lead can depend 
on the state of the system, |0) or |1). To illustrate the 
effect this dependence has on the cross correlations, we 
consider the case 


Majorana resonance F has to be smaller than the energy 
of the first subgap state. Importantly, the temperature 
T does not have to be smaller than F. 

In contrast, for the case of an accidental low-energy 
ABS, Prl is nonuniversal. In particular, it is sensitive to 
details such as the realization of disorder. 

The result of this work for the current cross correlation 
has its roots in the defining properties of MBSs. The 
high-voltage behavior can be shown to stem from the 
nonlocal nature of MBS; the fact that the occupation 
of the Majorana mode cannot be revealed by any local 
probe. The low-voltage behavior stems from the fact that 
the MBS induces perfect Andreev reflection at zero bias. 


ACKNOWLEDGEMENTS 


ro=0 ; r^ = 0 (19) 
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Appendix A: Hamiltonian Approach 


(VrVl) - (iVR)(iVL) = \ np { 1 - p ) = ^ (^1 
Inserting this into Eq. (15) one has 



le^~ 

~ 



( 21 ) 


The results presented in Eqs. ( 8 ) and (9) of Sec. HI can 
be derived from a Hamiltonian approach of transport. 
We start from an effective low-energy Hamiltonian de¬ 
scribing a multiple number of conducting channels which 
are coupled to a single MBS. Each of the channels belongs 
either to the left lead or to the right lead (although the 
calculation proceeds similarly in the case of a different 
number of leads). The Hamiltonian reads 


where (7 is a constant of order unity. Prl is monoton- 
ically increasing, asymptotically approaching a positive 
constat. This is in agreement with Eig. 3(b) and with 
the results of Ref. [35] . 


H^Hl + Ht, 

Hl = Y^= ii'^iKiUk + h . c .), 

ik ik 

(Al) 


VI. CONCLUSIONS 

When current from a topological superconductor is 
split into two metallic leads, the current cross correla¬ 
tion Prl exhibits universal behavior as a function of bias 
voltage V. The cross correlation is negative for all V and 
approaches zero at high voltage as —1/V. This behavior 
is robust and does not rely on a specific realization of the 
topological superconductor hosting the Majorana, or on a 
specific form of coupling to the leads. It can be observed 
even in disordered multichannel systems at finite tem¬ 
perature. For the effect to be observed the width of the 


where 7 describes the MBS, rjh, creates an electron with 
momentum k and energy eik in the i*** channel, and Xi is 
the coupling constant of the channel to the Majorana. 

In the wide-band limit the reflection matrix can be 
obtained by^^^a 


rtot(e) = 1 - (e + ittWmWI) ' Wm, (A2) 

with Wm being a vector of coupling constants given by 


{WM)i = \/Vi 


A. 

a: 


i = 1,...,4M 
i = 4M-kl,...,8M ’ 


(A3) 
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where Vi is the density of states of the i**' channel at the 
Fermi energy, and M is the number of spinful channels in 
each lead (all together there are 4M electronic channels). 
One obtains for the blocks of rtot [see also Eq. (2)] 


' i] ~ 


2TT^Xt\^ 
ie — T 


„he _ 
ij 


2TT^ViVj\i\j 


ie — T 


, (A4) 


with = [r®®(-e)]* and r®'‘(e) = [r^®(-e)]*, and 

where we have defined F = 2TTj2i=i 

Inserting Eq. (A4) into Eq. (7) results in 


dl 2 e 2 F^ 
dV ~ IT (eE)2 +r2’ 


(A5) 


and 




2e\ eV 

h ^ ^(eE) 2 +r 2 ’ 


(A 6 ) 


where F^ = We have therefore rederive 

Eqs. ( 8 ) and (9). We note that the definition of F here 
is in terms of the coupling constant, while in Sec. Ill it is 
given in terms of transmission amplitudes. In both cases, 
however, it equals the width of the Majorana-induced 
resonance. 


where Wl and IFr described the coupling to the left and 
right lead, respectively. As depicted in Fig. 5, each lead 
is coupled only to those lattice sites which are adjacent 
to it. Moreover, the coupling to each site is modulated 
according to the transverse profile of the particular chan- 



FIG. 5: Illustration of the tight-binding model corresponding 
to the system depicted in Fig. 1(a). Each lead is tunnel- 
coupled to the sites adjacent to it. The purple sites are ones 
in which there is a nonvanishing induced pair potential [cf. 
Eq. (11)]. 


nel. This is described by 


Appendix B: Details of Numerical Simulations 

To obtain the scattering matrix using Eqs. (11-13) we 
express the Hamiltonian H in first quantized form using 
a iN^Ny X AN^Ny matrix "Htb defined by 

H = Emn 'I'Utb'I'„ ; 4/t = ($t, $) , (Bl) 

where , ■, creates 

ra—2Ny {rix — l)+2(ny — l) + s r—{nxa,nya),s 

an electron with spin s on site (n^,ny) of an N^. x Ny 
square lattice. Here, s = 1 for spin =j' and s = 2 for spin 
j,. In our simulations we used = 90, and Ny = 6. 

The matrix W in Eq. (12) describes the coupling be¬ 
tween the extended modes of the leads and the sites of 
the lattice. In each lead there are M spinful transverse 
channels. In our simulations M = 4 (see Fig. 5). Includ¬ 
ing both leads, both spin species, and the particle-hole 
degree of freedom, IT is a AN^Ny x 8M matrix of the 
following form 

^ -w:) = > ^® 2 ) 


ITl = IT° (g) 



1 


0 

" (g) 0-0 . Wn = W°(g) 


0 

\o) 

Ny 

\y 




Wm sin 


0 


Trnm 

M+1 


1 < n < M 
M < n < Nx 


m = 1,..., M, 
(B3) 


where cr^ is a 2 x 2 identity matrix in spin space, and Wm is 
a set of coupling constants for each transverse channel of 
the leads. In this work we have used = 0.03Ao, Vm G 
{1,2,3,4}. 

Given the coupling matrix W and the first-quantized 
Hamiltonian Htb, the reflection matrix is calculated us¬ 
ing Eqs. (12) and (13). 
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